COUNTING ALL EQUILATERAL TRIANGLES IN {0,1, ...,n} 3 

EUGEN J. IONASCU 



Abstract. We describe a procedure of counting all equilateral triangles in the three dimensional 

space whose coordinates are allowed only in the set {0, 1, n}. This sequence is denoted here by 

ET(n) and it has the entry A102698 in "The On-Line Encyclopedia of Integer Sequences". The 

procedure is implemented in Maple and its main idea is based on the results in (3J. Using this we 

calculated the values ET(n) for n=1..55 which are included here. Some facts and conjectures about 

, ™ .„,.,,. In ET(n) 
this sequence are stated, ine mam of them is that hm exists. 



1. INTRODUCTION 

If we restrict the vertices of an equilateral triangle to be in Z 3 we obtain a typical element in 
£T(Z). It is not that hard to see that there are no such triangles whose vertices are contained in 
the coordinate planes or any other plane parallel to one of them. Also, the sides of a triangle in 
£T(Z) cannot be of an arbitrary length. If one such triangle is considered a whole family in £T(Z) 
can be generated from it that have vertices in the same plane. Moreover, we have shown in [3] the 
following theorems that we are going to use in our construction here. 

Theorem 1.1. If the triangle AOPQ G £T(Z) with O the origin and I = \ \OP\\ then: 

(i) the points P and Q are contained in a plane of equation ax + by + cz = 0, where 

(1) a 2 + b 2 + c 2 = 3d 2 , a, b, c, d £ Z 

and I 2 = 2q; 

(ii) the side length, I, is of the form y / 2(m 2 — mn + n 2 ) with m,n £ Z. 
It is important to be able to generate all the solutions of ([T|): 

Theorem 1.2. The following formulae give a three integer parameter solution of |7p: 

a = —x\ + x\ + x 2 — 2x\X2 — 2x1X3 
b = x \ - x\ + x\ - 2x x xi - 2x 2 x 3 

2,2 2 o n , Xl,X 2 ,X 3 G Z. 

C = X\ + Xg — X3 — ZX3X1 — ZX3X2 

d = x\ + x\ + x| 

Moreover, every solutions of (QJ), a, b, c, d can be found with simplifying if necessary by a 
common divisor of a, b, c and d with x\, x%, X3 no bigger in absolute value than \J ' 3 ~^ q. 
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We include some general remarks about the solutions of JT]) which are discussed in [3]: 

• if we assume that gcd(a, b, c) = 1 then all a, b, c, d must be odd integers 

• for every d odd there exist at least one solution which is not trivial (i.e. a = b = c = q) 

• [Gauss] positive integer n can be written as a sum of three squares iff n is not of the form 
4 k (8l + 7) with k, I € Z (see [1] for an elementary proof) 

Our construction depends essentially on a particular solution, (r,s) G Z 2 , of the equation: 



(3) 



2(a 2 + b 2 



s 2 + 3r 2 . 



It turns out that this Diophantine equation has always solutions if a, b, c and d are integers 
satisfying (pQ). The family we have mentioned can be described as another parametrization. 



Theorem 1.3. Let a, b, c, d be odd positive integers such thata 2 +b 2 



3d 2 , with gcd(d, c) 



1. Then a triangle AOPQ E £7~(Z) has its vertices in the plane of equation aa + b(3 + cry = iff 
P(u,v,w) and Q(x,y,z) are given by 



(4) 



with 



(5) 



u - 

V - 

w 



m u m — n u n, 
m v m — n v n, 
- m w m — n w n, 



and 



m x m — n x n, 



m z m — n z n, 



v 



— (rac + dbs)/q 
(das — bcr)/q 
r 



m x = — ^[db(3r + s) + ac(r — s)]/q, n a 
m y = i[cfa(3r + s) — bc(r — s)]/q, n 
m z = (r - s)/2, n z 
and 

m u = —(rac + dbs)/q, n u = — \[db(s — 3r) + ac(r + s)]/q 
m v = (das — rbc)/q, n v = ^[da(s — 3r) — bc(r + s)]/q 
m w =r, n w = (r + s)/2 



where q = a 2 + b 2 and (r, s) is a suitable solution of and m,n £ Z. 



Moreover, the side-lengths of this triangle are equal to dy / 2(m 2 — mn + n 2 



Let us observe that one can use Theorem 11.31 as long as min(gcd(q,c), gcd(q,a), gcd(d,b)) = 1. 
We have found the following counterexample to this property: a = 55063, b = 2396393, c = 
5(71)(2017)(1694953) and d = 3(41)(3361)(1694953) with gcd(a,d) = 41, gcd(b,d) = 3361 and 
gcd(c, d) = 1694953. However we have calculated that this property holds true for all solutions a, 
b, c, d of ([T]) such that gcd(a, b, c) = 1 and with all odd d < 4096. This allows us to calculate ET(n) 
for n = 1, .., 64 as we will see later. 
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2. Description of the procedure and some ingredients 

The idea is based on the facts above and a few other results. One would like to first find the 
side lengths of the triangles in £7~(Z) n C n . This will partition these triangles into clear classes. 
For this purpose we will use the Proposition 11.21 Then for a given side-length, I, we need to find 
all the possible planes that contain triangles of sides I . This gives another criteria of sub-partition 
even further these triangles. Using the parametrizations given in Theorem 11.31 then one finds the 
smallest such triangle within a given plane and which can fit in C n after a translation. Once that is 
obtained we have to rotate and translate it in all possible ways, but in a pairwise disjoint manner 
to fill out C n . A formula for the number of all these will be given. Finally all these numbers will 
be added up to make up ET(n). 

The first fact that we will use is the following geometric observation. 

Proposition 2.1. The largest side length of an equilateral whose vertices are contained in a 
cube of side lengths r is ry2. 

Proof. If an equilateral triangle of side lengths I, has its vertices in the cube [0, r] 3 , we denote 
by A\, A2, and ^.3 the areas of the three projections of this triangle on the three coordinate planes. 
It is easy to see that A 2 = A\ + A\ + A\ where A is the area of the equilateral triangle given. The 
maximum of the area of a triangle inscribed in a square of side lengths r is easy to see that is at 
most r 2 /2. Hence A 2 = 3/ 4 /16 < 3r 4 /4. This gives I < ry/2. Certainly this happens when the 
vertices are at the corners diagonally opposite on each face of the cube. ■ 
Let us work out a concrete example example: n = 4. Using Proposition 12.11 and the part (ii) 
of Theorem [TTJ the side lengths can be only y/2, y/6, 2y/2, y/U, 3^2, 2y/6, V26 or 4^2. The d 
values here are 1 or 3. Since 3(1 2 ) = l 2 + l 2 + l 2 and 3(3 2 ) = l 2 + l 2 + 5 2 are the only solutions 
of (H]), the parametrizations we need in this case are, as shown in [3]: 

^1,1,1 = {[(0, 0, 0), (to, —n, n — to), (to — n, —to, n)] : m, n G Z, to 2 + n 2 / 0} 

and 

^1,1,5 = {[(0, 0, 0), (4m — 3n, to + 3n, —to), (3m + n, —3m + An, —n)\ : to, n G Z, to 2 + n 2 ^ 0}, 

here we used the notation T a ^^ c standing for all triangles in £T(Z) having a vertex the origin and 
the other two in the plane {(a, j3, 7) : act + bf3 + 07 = 0}. 

Using the first parametrization we find the m, n such that the triangle obtained after a trans- 
lation fits in C4 and the side lengths are y/2: T\ = {(1, 0, 0), (0, 1, 0), (0, 0, 1)}. This triangle 
can be translated in various ways inside of C4, and together with all its cube symmetries and 
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translations contribute with a total of 512 in £T(4). We will prove a formula that gives the to- 
tal of all these triangles generated by T\ inside of C n . This parametrization has to be used for 
all the side lengths y/E, 2y/2, \/l4, 3\/2, 2\/6, \/26 and 4\/2: the corresponding triangles re- 
spectively are T 2 := {(1, 0, 2), (2, 1, 0), (0, 2, 1)}, 2T U T 3 := {(2, 0, 3), (0, 3, 2), (3, 2, 0)}, 3T 1; 2T 2 , 
T4 := {(1, 4, 0), (4, 0, 1), (0, 1, 4), }, 4T\. Using the same formula we will see the contribution of all 
these to ET(4) is respectively: 216, 216, 128, 64, 8, 16 and 8. 

There is need to use the second parametrization too since one can take d = 3 to obtain the side 
length 3\/2- This gives still a new triangle T 5 := {(0, 0, 1), (1, 4, 0), (4, 1, 0)} with a total of 96 other 
generated by it in C4. Tallying all these we get ET(A) = 1264. 

As we can see from this example, one has to derive a way of finding how many other triangles 
can a given one, say T, generate inside of C n under all possible translations and cube symmetries 
roughly speaking. We need to make this a little more precise. We are going to assume that the 
given triangle, T, that is inside Ct is minimal in the sense that at least one of the coordinates of 
the vertices in T is zero and t is the smallest dimension k of a cube containing T. 

Let us denote by 0(T) the orbit generated by T within C t under all translations and cube 
symmetries. We also need to introduce the standard unit vectors e\ = (1,0,0), e 2 = (0,1,0) and 
e 3 = (0,0,1). 

It is actually surprising that in order to compute the number /, of all distinct triangles generated 
by T (union of all translations of 0(T)) within C n in terms of the following five variables: 

(i) n -the dimension of the cube, 

(ii) t -the maximum of all the coordinates in T, 

(iii) a(T) -the cardinality of 0(T), 

(iv) f3{T) -the cardinality of 0{T) n [0{T) + ei], 

(v) 7 (T) -the cardinality of [0{T) + a] n [0(T) + e 2 ]. 

Theorem 2.2. The function f(T,n) described above is given by 
(6) f(T, n) = (n + 1 - tfa - 3(n + 1 - tf(n - t)(3 + 3(n + 1 - t)(n - tfj, 

for all n > t. 

Proof. Let us consider the cube C s = {0, 1, ...,s} 3 where s = n — t. Clearly the number of 
points in this set is (s + l) 3 . Each point p in the set C s is considered in here as a vector. So, 

/ = I U 0(T) +p\. One essential observation here is that \[0(T) +p\D [0(T) + q]\ = for every 

pec s 

p, q such that \\p — q\\ > 1, where \ \p — q\ \ = rain (\pi — qi\). This is due to the minimality of T. 

i=l,2,3 

Let us write the elements of C s in lexicographical order: pi,p 2 , ■■■■,Pk where k = (s + l) 3 . We look 
now at C s as the three dimensional grid graph. Faces in this graph are simply unit squares formed 
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with vertices from C s . One can look at the cardinality of U 0(T) + pi and show by induction on 

i=l..j 

j that this is equal to 

ja - (#edges(C s (j))(3 + (#faces(C s (j))j 

where C s (j) is the graph induced in C s by the vertices pi,P2, ■■■,Pj- Hence we just need to compute 
the number of edges and faces in C s . There are eight vertices in this graph that have degree 3 (the 
corners), there are (s — l) 3 vertices of degree 6, also 6(s — l) 2 vertices with degree 5 and finally 
12(s — 1) of degree 4. This gives a total of 

-[24 + 6(s - l) 3 + 30(s - l) 2 + 48(s - 1)] = 3s(s + l) 2 

edges. The number of faces is equal to ^[6s 3 + 6s 2 ] = 3s 2 (s + 1). ■ 
Example: Suppose we take T = T5. Then clearly t = 4. One can use a symbolic calculator to 

find a(T) = 96, (3(T) = 24 and j(T) = 0. So, the contribution of T5 to a cube C n if f(T^,n) = 

96(n - 3) 3 - 72(n - 3) 2 (n - 4) = 24n(n - 3) 2 . 

Remark: These facts give a way to find lower bounds for ET(n). For instance, if we put the 

contribution of T\ and T2 we obtain ET(n) > 8(2n — l)(n 2 — n + 1) for all n > 2. 

To generate the side lengths we would like use a well-known result due to Euler (see [1], pp. 568 

and [21, pp. 56). 

Proposition 2.3. [Euler's 6k + 1] An integer t can be written as m 2 — mn + n 2 for some 
m,n £ 7L if and only if in the prime factorization of t, 2 and the primes of the form 6k — 1 appear 
to an even exponent. 

3. The code 

Using Proposition 12.31 and Proposition 12 . 1 1 we have the following procedure in Maple to compute 
the side lengths modulo a factor of two and the square root: 

• sides:=proc(n) 

• local i,j,k,L,a,m,p,q,r,ms; 

• L := {l};ms := n 2 ; 

• for i from 2 to ms do 

• a := ifactors(i); k := nops(a[2]);r := 0; 

• for j from 1 to k do m := a[2][j][l]; p := m mod 6; q := a[2][j][2] mod 2; if r=0 and (m = 2 
or p = 5) and q = 1 then r := 1 fi; 

• if r=0 then L := L union {i};fi; 

• od; L := convert(L,list); end : 
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This procedure gives for n = 10: [1, 3, 4, 7, 9, 12, 13, 16, 19, 21, 25, 27, 28, 31, 36, 37, 39, 43, 
48, 49, 52, 57, 61, 63, 64, 67, 73, 75, 76, 79, 81, 84, 91, 93, 97, 100]. This gives the corresponding 
side-lengths 

[a/2, a/6,2a/2, a/II, 3^/2, 2\/6, a/26, 4\/2, a/38, a/42, 5a/2, 3a/6, 2\/l4, a/62, 6\/2, a/74, \/78, 
a/86, 4\/6, 7\/2, 2\/26, \/Tl4, \/l22, 3\/l4, 8\/2, a/134, \/i~46, 5\/6, 2\/38, \/l58, 9\/2, 2\/42, \/l82, 
\/l86, a/194, 10\/2] 

We need a procedure that will give the odd values of d that "divide" a certain side length in the 
sense it is possible to write it as dV m 2 — run + n 2 with m,neZ: 

• cZfcZ := proc(side) 

• local i,x,noft,div,y,yl,z; 

• x := convert(divisors(side),list);noft := nops(x);div := {}; 

• for i from 1 to noft do z := mod 2; 

• if z = 1 then y := side/x[i] 2 ; yl := floor(y); if y = yl then div := div union h;h; 

• od; 

• convert (div, list); end: 

For instance, if side = a/882 this procedure gives [1,3,7,21] which means we have at least four 
possible parametrizations that we can use to find minimal equilateral triangles. 

Next we need to find all the nontrivial solutions [a, b, c] of (pQ), given and odd positive integer d, 
with the property gcd(a, b, c) = l, < a < b < c which is based on an internal procedure to solve 
Diophantine equation A = X 2 + Y 2 : 

• abcsol:=proc(q) local i,j,k,u,x,y,sol,cd; sol := {}; 

• for i from 1 to q do 

• u := [isolve{2>d 2 — i 2 = x 2 + y 2 )]; k := nops(u); 

• for j from 1 to k do 

• if r/is(u[j][l]) >= % and rhs{u[j] [2]) >= i then 

• cd := gcd(gcd(i,rhs(u[j][l])),rhs(u[j][2])); 

• if cd = 1 then sol := sol union {sort([i, rhs(u[j] [1]), rhs(u[j] [2])])};fi;fi; 

• od;od; convert (sol, list); end: 

For d = 17, ancsol finds four different solutions, [[11, 11, 25], [13, 13, 23], [1, 5, 29], [7, 17, 23]], and 
in a few seconds sends out 333 solutions for d = 2007. One interesting solution in this last case is 

1937 2 + 1973 2 + 2107 2 = 3(2007) 2 . 

Now based on the Theorem 11.31 we take a solution of ([T]) as given by the procedure above and 
calculate the general parametrization: 
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• findpar:=proc(a,b,c,m,n) 

• local ij,r,s,sol,mx,nx,my,ny,mu,nu,mv,nv,mz,nz,mw,nw,q,d,u,v,w,x,y,z,ef,ns,om,l; 



• q := a 2 + b 2 ; sol := convert(isolve(2q = x 2 + 3y 2 ),list); ns := nops(sol); d := y « 2 +fe 2 +c 2 ■ 
ef := 0; 

• for i from 1 to ns do 

• if ef = then r := r/is(soZ[«][l]); s := rhs(sol[i][2]); 

• if s 2 + 3r 2 = 2g then mz := (r — s)/2; nz := r; mw := r; nw := (r + s)/2; rax : = 
— (db(3r + s) + ac(r — s))/(2q); nx := — (rac + dbs)/q; my := (da(3r + s) — 6c(r — s))/(2q); 
ny := —{rbc — das)/q; mu := —(rac + dbs)/q; nu := —(db(s — 3r) + ac(r + s))/(2q)\ 
mv := (das — rbc)/q; nv := —(da(3r — s) + 6c(r + s))/(2q); 

• if ira = floor (mx) and nx = floor(nx) and my = floor {my) and ny = floor (ny) and 

= floor(mu) and nu = floor (nu) and mi) = floor (mv) and = floor (nv) then 

• n := (mu)m — (nu)n; v := (mv)m — (nv)n; w := (mw)m — (nw)n; x := (mx)m — (nx)n; 
y := (my)m — (ny)n; z := (mz)m — (nz)n; om := [[u, v, w], [x, y, z}]; 

• ef := 1; fi;fi; fi; od; om; end: 

For the solution, [1, 5, 29], found earlier for the case d = 17, findpar gives 

[[-11m - 13n, -21m + 20n, 4m - 3n], [-24m + lln, -m + 21n, m - 4n]]. 

Next, using this parametrization we would like to find if there is any equilateral triangle in £7~(Z) 
which after a translation fits inside C stopp . 

1: minimaltr:=proc(s,&,b,c,stopp) 

2: local i,z,u,nt,d,m,n,T,a,/?,7,tr,out,L,tri,noft,tria,orb,avb, length, lengthn; 

3: d := y/ a2 +fr 2 + ii; z := ^.j u := convert(isolve(z = q 2 — qr + r 2 ),list); nt := nops(u); 

4: for i from 1 to nt do 

5: T := findpar(a, b, c, rhs(u[i][l]),rhs(u[i][2])); 

6: a:=mm(T[l][l],T[2][l],0); /3 := mm(T[l][2],T[2][2],0); 7 == mm(T[l][3],T[2][3],0); 

7: tr[i] := [T[l][l] - a, T[l][2] - /?, T[l][3] - 7 ], [T[2][l] - a, T[2][2] - /3, T[2][3] - 7 ], [-a, - 

out[i] := max(tr[i] [1] [1], tr[i\ [1] [2], tr[i] [1] [3], tr[i\ [2] [1], 
8 - trH[2][2],trH[2][3],trH[3][l],trH[3][2],tr[z][3][3]); 
9: od; L := sort([seg(oui[i], i = l..nt)]); tri := {}; 

10: for i from 1 to nt do if out[i] <= stopp then tri := tri union {ir[i]}; fi; od; 
11: tri := convert(tri,list); tria := {}; 

12: if nops(tri) > then noft := nops(tri); tria := |in[l]}; orb := trans^(tri[l]); 
13: for i from 1 to noft do avb := evalb(tri[i] in orb); 

14: if at>6 = false then or6 := or& union transl(tri[i]); tria := tria union {irz[z]}; 
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15: fi; od; fi; tria; end: 

The minimal triangle given by this procedure for s = 17s/2, a = 1, b = 5, c = 29, stopp = 30: 
{[11, 21, 0], [24, 1, 3], [0, 0, 4]}. The last part of the procedure is actually searching for a set of 
triangles that generate all the triangles in £7~(Z) that lay in planes of normal (a, b, c) or all other 
23 possibilities obtained by permuting the coordinates and changing signs. The next procedure is 
used above and later in order to compute the parameters a(T), @{T) and 7(T). 

• transl:=proc(T) 

• local S,Q,i,j,k,a2,b2,c2,a,b,c,d; 

• Q :=convert(T,list); a := max(Q[l][l], Q[2][l], Q[3][l]); b := max(Q[l}[2], Q[2][2], Q[3][2]); 
c := max(Q[l] [3], Q[2] [3], Q[3] [3]); d := max(a, b, c); a2 := d — a; b2 := d — b; c2 := d — c; 
S ■= orbit(T); 

• for i from to a2 do 

• for j from to b2 do 

• for k from to c2 do 

• S := S union orbit(addvect(T, [i,j, k])); 

• od;od;od;S;end: 

Here the procedure addvect and orbit are: 

1: addvect :=proc(T,v) local Q,a,b,c; 

2: Q := convert(T,list); a := v[l}; b := v[2]; c := v[3]; 

3: { [Q [1] [1] + a, Q[l] [2] + b,Q[l] [3] +c], [Q [2] [1] +a,Q [2] [2] +b,Q [2] [3] +c],[Q [3] [1] +a,Q [3] [2] + 

6,Q[3][3]+c]}; 
4: end: 

and 

• orftit := proc(T) 

• localS, Q,T1; 

• Q := convert(T, list); 

. Tl :={[Q[l][3],Q[l][2],g[l][l]],[Q[2][3],Q[2][2],Q[2][l]],[Q[3][3],Q[3][2],Q[3][l]]}; 

• 5 : = orbit\{T) union or6itl(Tl); 5; 

• end: 

The orbitl takes care of the cube symmetries: 

• orbitl:=proc(T) local 

• i,k,Tl,a,b,c,x,T2,T3,T4,T5,T6,T7,T8,T9,T10,Tll,T12,T13,T14,T15,T16,T17,T18, 

• T19,T20,T21,T22,T23,T24,S,Q,d,al,bl,cl; Q:=convert(T,list); 
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a: 



max(Q[l][l],Q[2][l],Q[3][l]); a l:=min(Q[l][l],Q[2][l],Q[3][l]); 
b:=max(Q[l] [2],Q[2] [2],Q[3] [2]);bl:=min(Q[l] [2],Q[2] [2] ,Q[3] [2]); 
c:=max(Q[l][3],Q[2][3],Q[3][3]);cl:=min(Q[l][3],Q[2][3],Q[3][3]); 
d:=max(a,b,c); T1:=T; 

T2:=[Q[1][2],Q[1][3],Q[1][1]],[Q[2][2],Q[2][3],Q[2][1]],[Q[3][2],Q[3][3],Q[3][1]]; 
[Q[1][1],Q[1][3],Q[1][2]],[Q[2][1],Q[2][3],Q[2][2]],[Q[3][1],Q[3][3],Q[3][2]] ; 
[Q[l][l],Q[l][2],d-Q[l][3]],[Q[2][l],Q[2][2],d-Q[2][3]],[Q[3][l],Q[3][2],d-Q[3][3]] 



T3 
T4 
T5 
T6 
T7 
T8 
T9 
T10 
Til 
T12 
T13 
T14 
T15 
T16 
T17 
T18 
T19 
T20 
T21 
T22 
T23 



= [Q[l][2],Q[l][3] ) d-Q[l][l]],[Q[2][2] 
:[Q[l][l],Q[l][3],d-Q[l][2]],[Q[2][l] 
=[Q[l][l],d-Q[l][2],Q[l][3]],[Q[2][l] 
=[Q[l][2],d-Q[l][3],Q[l][l]],[Q[2][2] 
= [Q[l][l],d-Q[l][3],Q[l][2]],[Q[2][l] 



=[d-Q[l][l],Q[l] 
=[d-Q[l][2],Q[l] 
=[d-Q[l][l],Q[l] 

=[Q[l][l],d-Q[l] 
=[Q[l][2],d-Q[l] 

=[Q[l][l],d-Q[l] 
=[d-Q[l][l],d-Q 

=[d-Q[l][2],d-Q 
=[d-Q[l][l],d-Q 
=[d-Q[l][l],Q[l] 
=[d-Q[l][2],Q[l] 
=[d-Q[l][l],Q[l] 

=[d-Q[l][l],d-Q 
=[d-Q[l][2],d-Q 
=[d-Q[l][l],d-Q 



[2],Q[l][3]],[d-Q[2 
3],Q[l][l]],[d-Q[2 
3],Q[l][2]],[d-Q[2 



d-Q[l][3] 
d-Q[l][l] 
d-Q[l][2] 
2],Q[1][3] 

3],Q[1][1] 
3],Q[1][2] 
d-Q[l][3] 
d-Q[l][l] 
d-Q[l][2] 
2],d-Q[l] 

3],d-Q[l] 
3],d-Q[l] 



,[Q[2 
,[Q[2 
,[Q[2 



Q[2][3],d-Q[2][l]] 
Q[2][3],d-Q[2][2]] 
d-Q[2][2],Q[2][3]] 
d-Q[2][3],Q[2][l]] 
d-Q[2][3],Q[2][2]] 
[1],Q[2][2],Q[2][3 
[2],Q[2][3],Q[2][1 
[1],Q[2][3],Q[2][2 
[l],d-Q[2][2],d-Q 
[2],d-Q[2][3],d-Q 
[l],d-Q[2][3],d-Q 
,[d-Q[2][l],d-Q[2][2],Q 
,[d-Q[2][2],d-Q[2][3],Q 
,[d-Q[2][l],d-Q[2][3],Q 
,[d-Q[2][l],Q[2][2],d-Q 
,[d-Q[2][2],Q[2][3],d-Q 
,[d-Q[2][l],Q[2][3],d-Q 



[Q[3][2],Q[3][3],d-Q[3][l]] 

[Q[3][l],Q[3][3] s d-Q[3][2]] 

[Q[3][l],d-Q[3][2],Q[3][3]] 

[Q[3][2],d-Q[3][3],Q[3][1]] 

[Q[3][l],d-Q[3][3],Q[3][2]] 

|,[d-Q[3][l],Q[3][2],Q[3][3]] 

|,[d-Q[3][2],Q[3][3],Q[3][l]] 

|,[d-Q[3][l],Q[3][3],Q[3][2]] 



2] [3] 
2][1] 
2] [2] 
2] [3] 

2][1] 
2] [2] 
2] [3] 

2][1] 
2] [2] 



,[Q[3][l],d-Q[3][2],d-Q[3][3] 
,[Q[3][2],d-Q[3][3],d-Q[3][l] 
,[Q[3][l],d-Q[3][3],d-Q[3][2] 
,[d-Q[3][l],d-Q[3][2],Q[3][3] 
,[d-Q[3][2] s d-Q[3][3],Q[3][l] 
,[d-Q[3][l] s d-Q[3][3],Q[3][2] 
,[d-Q[3][l],Q[3][2],d-Q[3][3] 
,[d-Q[3][2],Q[3][3],d-Q[3][l] 



[d-Q[3][l],Q[3][3],d-Q[3][2] 
3]],[d-Q[2][l],d-Q[2][2],d-Q[2][3]] ) [d-Q[3][l],d-Q[3][2],d-Q[3][3]] 
l]],[d-Q[2][2],d-Q[2][3],d-Q[2][l]], [d-Q[3][2],d-Q[3][3],d-Q[3][l]] 
2]],[d-Q[2][l],d-Q[2][3],d-Q[2][2]], [d-Q[3][l],d-Q[3][3],d-Q[3][2]] 



T24 

S := {Tl, T2, T3, T4, T5, T6, T7, T8, T9, T1Q, Til, T12, T13, T14, T15, T16, T17, T18, T19, T20, 
T21, T22, T23, T24}; 
S; end: 



Finally, we are ready to calculate the parameters in Theorem 
(3(T) = inters(T) where 



We have ct(T) = transl(T), 



inters:=proc(T) local a,b,c,Q,d,S,m,i,Sl,S2; Q:=convert(T,list); 
a:=max(Q[l][l],Q[2][l],Q[3][l]); 
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. b:=max(Q[l][2],Q[2][2],Q[3][2]); 
. c:=max(Q[l][3],Q[2][3],Q[3][3]); 

• d:=max(a,b,c);S2:=transl(T);S:=convert(S2,list);m:=nops(S);Sl:=; 

• for i from 1 to m do 51 := 51 union {addvect(S[i], [0,0, 1])}; od; 

• 52 intersect 51; end: 
and j(T) = intersch(T) where 

• intersch:=proc(T) local a,b,c,Q,d,S,m,i,Sl,S2,S3,S4; 

• Q:=convert(T,list); 

• S2:=transl(T);S:=convert(S2,list);m:=nops(S);Sl: = ; 

• for i from 1 to m do 

• 51 := 51 union {addvect(S[i], [0, 0, 1])}; od; 53 := {}; 

• for i from 1 m do 53 := 53 union {addvect(S[i], [0, 1,0])}; od; 

• nops(Sl intersect S3); end: 

The Theorem 12.21 is then implemented in 

/ := (n, d, a, (3, 7) -> (n - d + lfa -3(n-d+ l) 2 (n - d)(3 + 3(n - d + l)(n - dfj 

• notrincn:=proc(T,n) 

• local Q,a,b,c,x,a2,b2,c2,d,y,z,w; 

• Q:=convert(T,list); 

. a2:=max(Q[l][l],Q[2][l],Q[3][l]); 
. b2:=max(Q[l][2],Q[2][2],Q[3][2]); 
. c2:=max(Q[l][3],Q[2][3],Q[3][3]); 

• d:=max(a2,b2,c2); 

• x:=nops(transl(T));y:=nops(inters(T));w:=intersch(T); 

• z := f(n,d,x,y,w); end: 

In the end one has to put together all these procedures and add the number of triangles together. 

• main:=proc(p,lastside,nuptols) 

• local i,j,k,s,nos,div,nod,nop,sol,x,netr,noft,l,z; netr:=nuptols; 

• s:=sides(p);nos:=nops(s);print(s); 

• for i from lastside to nos do 

• div:=dkl(s[i]);nod:=nops(div); 

• for j from 1 to nod do 

• sol:=abcsol(div[j]);nop:=nops(sol); 

• for k from 1 to nop do 

• x:=minimaltr (s[i] ,sol[k] [1] ,sol[k] [2] ,sol[k] [3] ,p) ; 
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• noft:=nops(x); if no ft >= 1 then 

• for 1 from 1 to noft do 

• z := notrincn(x[l],p); 

• netr:=netr+z; print (s[i],div,sol[k],x[l],checkeq(x[l]),z,netr,i); 

• od; fi; od; od; od; netr; end: 

The values ET{n) for n = 1...55 computed with main are in given below in increasing order: 
8, 80, 368, 1264, 3448, 7792, 16176, 30696, 54216, 90104, 143576, 220328, 326680, 471232, 664648, 
916344, 1241856, 1655208, 2172584, 2812664, 3598664, 4553800, 5702776, 7075264, 8705088, 10628928, 
12880056, 15496616, 18523472, 22003808, 26000584, 30567400, 35756776, 41631672, 48278136, 
55753272, 64134536, 73495760, 83924408, 95513248, 108379264, 122661840, 138315720, 155613408, 
174622488, 195478424, 218279240, 243170376, 270288064, 299790968, 331832248, 366610560, 404253120, 
444911712, and 488902856. 

4. Some facts and conjectures 

If we look at the sequence a n = -j^+jj; n £ N it seems like it is increasing. This sequence is 
clearly bounded from above since the number of all triangles in {0, n} 3 is not more than (n + 1) 9 
and so a n < 9. Numerically, the best upper-bound for a n seems to be 5 which is equivalent to 
saying that ET(n) < (n + l) 5 for all n £ N. 

From what we have seen before each class of triangles determined by a, b, c, d, a solutions of ((H), 
brings in a contribution that is a polynomial in terms of n. If we add these polynomials together, 
we get a polynomial which can be expressed is in the variable C, = n — 1 (n = £ + 1) as follows: 

n = 1: pi(C) = 8C 3 + 24C 2 + 24C + 8, ET(l) = Pl (0) = 8; 

n = 2: p 2 (C) = pi(C) + 16 + 48(C - 1) + 16(C - l) 3 + 48(C - l) 2 , ET{2) = p 2 (l) = 80; 
n = 3: pg(C) = p 2 (C) + 24 + 72(C - 2) + 24(C - 2) 3 + 72(C - 2) 2 , ET(S) = p 3 (2) = 368; 
n = 4: p 4 (C) = Ps(0 + 128 + 312(C - 3) + 56(C - 3) 3 + 240(C - 3) 2 , ET(4) = p 4 (3) = 1264; 
n = 5: p 5 (C) = P4(C) + 4 + 120(C - 4) + 40(C - 4) 3 + 120(C - 4) 2 , ET{5) = p 5 (4) = 3448; 
n = 6: p6(0 = p 5 (0 + 4 8 + 144(C - 5) + 48(C - 5) 3 + 144(C - 5) 2 , ET{6) = Pe {5) = 7792; 
n = 7: p 7 (C) = Pe(C) + 776 + 1392(C - 6) + 128(C - 6) 3 + 744(C - 6) 2 , ET(7) = p 7 {6) = 16176; 
n = 8: p 8 (C) = p 7 (0 + 232 + 552(C - 7) + 88(C - 7) 3 + 408(C - 7) 2 , ET(8) = p 8 (7) = 30696; 
n = 9: pq(C) = p 8 (C) + 360 + 840(C - 8) + 120(C - 8) 3 + 600(C - 8) 2 , ET(9) = p 9 (8) = 54216; 
n= 10: pio(C) =p 9 (C) + 80 + 80(C-9) 3 + 240(C-9) 2 + 240(C-9), ET (10) = p w (9) = 90104; 



We conjecture that in general 
(7) Pn(C) = Pn-l(C) + «n(C - n + l) 3 + u n (C - n + l) 2 + w„(C - n + 1) + s n , n G N, 
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with u n ,v n ,w n , and s n non-negative integers. 



5 i 




10 20 30 40 50 



As the graph above of n — > ln( - w+1 ) suggests, the second conjecture is that the following limit 
exists 



(8) lim ^M = c . 

- m(n + 1) 



n—>oo 
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